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Abstract 

In the framework of the Littlest Higgs Model with T-parity (LHT), we study the contributions 
of the new particles to Zbb couplings at one-loop level. Based on these results, we further study 
the branching ratio Rh and the unpolarized forward- backward asymmetry A^p^. We find that 
the correction of the new particles to Zbb couplings is mainly on the left-handed coupling and 
has small part of the parameter space to alleviate the deviation between theoretical predictions 
and experimental values. The precision measurement value of Rb can give severe constraints on 
the relevant parameters. The constraints from the precision measurement value of A^p^ are very 
weak. 

PACS numbers: 14.65.Fy,12.60.-i,12.15.Mm,13.85.Lg 
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I. INTRODUCTION 



The Standard Model (SM) has been very successful, however, it is still believed to be a 
theory effective at the electroweak scale and some new physics (NP) must exist at higher 
energy regimes. So far there have been many speculations on the possible forms of the NP 
beyond the SM, one of the interesting possibilities is the Little Higgs model. The little 
Higgs theory was proposed [1] as a possible solution to the hierarchy problem and remains 
a popular candidate for the NP. The Littlest Higgs (LH) model [2i] is a cute economical 
implementation of the little Higgs, but suffered from severe constraints from electroweak 
precision tests j^, which would require raising the mass scale of the new particles to far 
above TeV scale and thus reintroduce the fine-tuning in the Higgs potential ^. The most 
serious constraints resulted from the tree-level corrections to precision electroweak observ- 
ables due to the exchanges of the additional heavy gauge bosons present in the theories, 
as well as from the small but non- vanishing vacuum expectation value (VEV) of an addi- 



tional weak-triplet sca^ 
T-parity is proposed 



ar field. In order to solve this problem, a discrete symmetry called 
^, which explicitly forbids any tree- level contributions from the 
heavy gauge bosons to the observables involving only the SM particles as external states. 
The interactions that induce triplet VEV contributions is also forbidden. This model is 
called the Littlest Higgs Model with T-parity (LHT). In the LHT model, corrections to 
the precision electroweak observables are generated exclusively at loop level. 

The branching ratio Rb is very sensitive to the NP beyond the SM, the precision 
experimental value of Rb may give a severe constraint on the NP gI . Experimentally, the 
electroweak observables have been precisely measured at the SLC and LEP, in the most 
recent analysis of the electroweak data, Rb = 0.21629 ±0.00066 differs from the SM fit by 
0.7(T, A^pjj = 0.0992 ± 0.0016 disagrees with the SM fit by -2.9cj Q- Furthermore, the 
experimental value of Zbb couplings disagrees with the SM fit by about 3a, especially the 
deviation of the right-handed coupling is so large that it is very difficult to explain. These 
significant deviations from the Ap^ and the Zbb couplings might be the first window 
into the NP. With the running of the LHC, they will be further researched. In the LHT 
model, there are new fermions and new gauge bosons, which can contribute to the Zbb 
couplings and give modifications to the Rb and A'^p^. Therefore, it is possible to give some 
constraints on the relevant parameters via their radiative corrections to the Rb and A''p^. 



In this paper, we calculate the contributions of the LHT model to the Zbb couplings. On 
this basis, we further study the Rb and A^p^, then we give the constraints on the relevant 
parameters according to the precision measurements. 

This paper is organized as follows. In Sec. II we recapitulate the LHT model and discuss 
the new flavor interactions which will contribute to the Zbb vertex. In Sec. Ill we calculate 
the one- loop contributions of the LHT model to the Zbb vertex, Rh and A^p^, then the 
relevant numerical results are shown. Finally, we give our conclusions in Sec. IV. 



II. A BRIEF REVIEW OF THE LHT MODEL 



The LHT [5[ is based on a non-linear sigma model describing the spontaneous breaking 
of a global SU{5) down to a global 5*0(5). This symmetry breaking takes place at the 
scale / ~ 0{TeV) and originates from the VEV of an SU{5) symmetric tensor S, given 
by 

/ \ 

02x2 l2x2 
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From the SU{5)/ S0{5) breaking, there arise 14 Goldstone bosons which are described by 
the "pion" matrix H, given explicitly by 
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Under T-parity the SM Higgs doublet, H = {-i7i+y/2, {v + h + i7r^)/2f is T-even while 
other fields are T-odd. 

The Goldstone bosons u^, u^, rj are respectively eaten by the new T-odd gauge bosons 
W^, Zh, Ah, which obtain masses at 0{v'^ / p) 



Mwn = Mz„ = gf{l 



V 



with g and g' being the SM SU{2) and f/(l) gauge couplings, respectively. 



(3) 



The Goldstone bosons 71^,7?*^ are eaten by the T-even Zi^ bosons of the SM, 

which obtain masses at 0{v'^ / p) 

Mw, = ^(1 - T^), Mz, = ^^(1 - T^) (4) 
2 12/^ 2 cos Ow 12/^ 

The photon Al is also T-even and remains massless. 

For each SM fermion, a copy of mirror fermion with T-odd quantum number is added 
in order to preserve the T-parity. For the mirror quarks, we denote them by m^, d^, where 
i= 1, 2, 3 are the generation index. At the order of 0{v'^ / p) their masses are given by 

,2 



y/2Kif, m^. = m^. (1 - ^) (5) 



m 



where the diagonahzed Yukawa couphngs of the mirror quarks. 

In order to cancel the quadratic divergence of the Higgs mass induced by top loops, 
an additional heavy quark is introduced, which is even under T-parity. The imple- 
mentation of T-parity then requires also a T-odd partner T~ . Their masses are given 
by 

where xl is the mixing parameter between the SM top-quark t and the new top-quark 
T+. 

Just like the SM, the mirror sector in the LHT model also has weak mixing, parame- 
terised by unitary mixing matrices: two for mirror quarks and two for mirror leptons: 

VHu,VHd,VHl,VH. (8) 

Vhu and Vhci are for the mirror quarks which are present in our analysis. Vhu and Vnd 
satisfy the physical constraints V^J/nd = ycKM- We follow js] to parameterize Vnd with 
three angles 9^2, ^23 ' ^13 three phases 5^2, ^2^, Sf^ 

'-12'-13 '^12'-13^ ■^IS^ 

""512^23^*'^''^ ~ C^2"523"5l3'2*^'^''^~'^^^'' ^^2^23 ~ sf2S23S^3e**^'^i3~'^i2~^23) SgsCfgC'^'^as 
y 5^25236*^^12+^23) — 0^26235^36*^13 — cf2523e*''23 — S^2C23'5l3e*^^i3~^i2) 6336^3 J 



V, 



Hd 
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III. THE ONE-LOOP CORRECTIONS TO Zhb COUPLINGS IN THE LHT 
MODEL 



We employ the following notation for the effective Zhh interaction: 

ow^W 

= T^^^mg'v-gh^bz,) (10) 

where 6w is the Weinberg angle, Sw = sin^vK, Cw = cosOw, Pl = and Pr = 
The effective couplings are then written as 

Ar = + Sgt% + 5gZ (n) 

Aa = Aa + ^9vA + (12) 

where gLR,gvA respectively the radiatively-corrected effective couplings, g'l are 
respectively the left-handed and right-handed Zbb couplings at tree level, Sgf'^ and Sg^^ 
are their corresponding one-loop corrections of the SM and the NP, gy^ are respectively 
the vector and axial vector coupling coefficients of Zbb interaction at tree level, Sgy^J^ and 
Sgy^^ are their corresponding one-loop corrections of the SM and the NP. The tree-level 
couplings are given by 

gv = g''L + gR = -l + ls'w , g\ = gi-g'k = -l (i4) 

The branching ratio is defined as 

^ TjZ ^ bb) 
T{Z ^ hadrons) 

The full hadron width is the sum of widths of five quark channels: 
T{Z -» hadrons) = T{Z -» uu) + T{Z -> dd) + T{Z ss) + T{Z cc) + T{Z bb) (16) 
For the decays to any of the five pairs of quarks qq we have[ 

r, = r(z ^ qq) = UToiglgRAg + gl.Rv,) (17) 
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with To = 24^27^ ' ■'^^^^ ^^"^ fi'v'q are the axial- vector and effective vector couphngs. 
The radiators Rxq and Ry^ contain contributions from the final state gluons and photons. 
In the crudest approximation 



R 



Vq 



Rao = 1 + 



a. 



where (y.s{(f) is the QCD running couphng constant: 



a. = a.. 



Mi 



(19) 



The expression of the radiative correction to Rb can be expressed as [10 1 



5Rb 



■[gvb{'i~P^)5gvh + '2gAbP^5g_ 
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with [5 



1 



4m: 



7^ being the velocity of b-quark in Z decay, here rhf, is the value of the 



running mass of the b-quark at scale Mz^ calculated in MS scheme 

The unpolarized forward-backward asymmetry in the decay to hb equals: 

Nf-Nb 



A 



FB 



(21) 



Nf + Nb 

where Nf is the cross section for finding the scattered fermion in the hemisphere defined 
by the incident electron direction and Nb is the cross section for finding it in the positron 
hemisphere. It can be expressed as 



-d 



TT 



(22) 



where the factor (1 — ^) represents a QCD radiative correction, as in Ref. [12|], for which 
we use the numerical value 0.95, refers to the creation of Z boson in e^e~ -annihilation, 
while Aft is the left-right coupling constant asymmetryrefers to its decay in bb [9|] 

'^gAhgvb 



Ah 



(23) 



P^gl, + (3 - f3^)glJ2 

The relevant Feynman diagrams for the LHT contributions are shown in Fig.l. We use 
the 't Hooft-Feynman gauge, so the contributions of Goldstone bosons should be involved. 
In our calculation, gAb and gyh should be replaced by gAb and gvb, S'yA + ^9v^a "^an be 



found in Ref. 



13| . The calculations of the loop diagrams are straightforward. Each loop 



diagram is composed of some scalar loop functions [14], which are calculated by using 



(d) 
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FIG. 1: Feynman diagrams of Z — t- 66 at one-loop level in the LHT model. 



LOOPTOOLS [15]. The relevant Feynman rules can be found in Ref. [16|. We applied 
the on-shell renormalization scheme and have checked that the divergences are canceled. 

In the numerical calculations we take the input parameters 17|] as Fermi constant 
Gp = 1.16637 X lO^^GeV^"^, the fine-structure constant a = 1/128, Z-boson mass Mz^ = 
91.26*6^^, fermion masses rrif, the electroweak mixing angle = 0.231 and the final-state 
asymmetry parameter Af, = 0.1515. In our calculation, the relevant LHT parameters are 
the scale /, the mixing parameter xl, the mirror quark masses and parameters in the 
matrices Vhu and Vnd- 

For the mirror quark masses, from Eq.(5) we get m^i^ = m^i^ at 0{v/f) and further 
assume 

Mi2,m„3 =mrf3 =M3 (24) 



For the matrices Vhu and Vnd, considering the constraints in Ref.[l8|, we study the 
completely generic scenario, i.e. the six parameters of Vnd are arbitrary. After that, we 



follow Ref.[l9| to consider the following two scenarios for comparison: 
Scenario I: V^d = 1, Vhu = Vckm 

6: 



Scenario II: S'}^ = 0.5, Sf^ 



-'23 



0, 



xSM qd 
'^13 ' '^ij 



5',^*^otherwise 



Firstly, we discuss the Rb changes with the LHT parameters, the numerical results are 
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FIG. 2: Scatter plots of Rb versus xiin arbitrary scenario, scenario I and scenario II, respectively. 
The experimental value = 0.21629 ± 0.00066. 
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FIG. 3: Scatter plots of Rh versus M12 in three different scenarios, respectively. The experimental 
value Rb = 0.21629 ± 0.00066. 

summarized in Fig. (2-3). To see the influence of the mixing parameter xl on the Rb, 
considering the existing constraints, we let the parameters vary randomly in the range: 
M12 = 300 ~ 3000GeF, M3 = 300 ~ 3000Gey, / = 400 ~ 3000Ge\/. In these three 
scenarios, we can see the plots of Rb decline with the xl increasing, which shows that the 
contribution of the mixing diagrams between t and T"*" is negative and becomes larger 
with the Xl increasing. When xl > 0.7, part of the plots are beyond the 2a regions of its 
experimental value. This feature is similar in three different scenarios. 

To see the influence of the first two generation mirror quarks mass M12 on the Rb, 
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FIG, 4: Scatter plots of A^pj^ versus xl in three different scenarios, respectively. The experi- 
mental value ^5^5 = 0.0992 ± 0.0016. 
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FIG. 5: Scatter plots of A^j,^ versus / in three different scenarios, respectively. The experimental 
value A'^pj^ = 0.0992 ± 0.0016. 



considering the constraint from Rb on the X/,, we let the parameters vary randomly in 
the range: M3 = 300 ~ SOOOGeV, / = 400 ~ SOOOGeV", xl = 0.1 ~ 0.7. In these three 
scenarios, we can see the plots of Rt are almost in the 2a regions of its experimental value. 
The noticeable feature is that the Rb isn't sensitive to M12 so that the constraint from Rb 
on M12 is very loose. 

Secondly, we discuss the A^p^ changes with the LHT parameters, the numerical results 
are summarized in Fig. (4-6). Same as the Rb, the A'^p^ isn't sensitive to M12, so we 
don't give the figures of the A'^p^ as the function of M12. To see the influence of the 
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FIG. 6: Scatter plots of A'^pj^ versus M3 in three different scenarios, respectively. The experi- 



mental value A},^ = 0.0992 ± 0.0016. 




FIG. 7: The left-handed and right-handed coupling constants in the LHT model. The experi- 
mental value g'l = -0.4182 ± 0.0015, 5^ = 0.0962 ± 0.0063. 




FIG. 8: The effective vector and axial-vector coupling constants in the LHT model. The exper- 
imental value glr = -0.3220 ± 0.0077,5^1 = -0.5144 ± 0.0051. 
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mixing parameter xl on the A^pBi parameters vary randomly in the range: 

Mi2 = 300 ~ 3000GeF, M3 = 300 ~ 3000Ge\/, / = 400 ~ 3000Gey. For the same 
reason, we can see the plots of A^p^ decline and become closer to the experimental central 
value with the xl increasing. However, the contribution of the new particles is not large 
enough so that the plots of the A^p^ are still entirely scattered between the 2a and 3cr 
region of its experimental value. 

To see the influence of the scale / on the A^^, we let the parameters vary randomly 
in the range: M12 = 300 ~ 3000GeV", M3 = 300 ~ 3000GeV", = 0.1 ~ 0.7. We can 
see the plots of A^p^ are entirely between the 2(T and 3(T region of its experimental value. 
The plots of A^pB become closer to the SM with the / increasing, which shows that the 
contribution of the heavy particles decouples with the / increasing. 

To see the influence of the third generation mirror quarks mass M3 on the A^p^, we let 
the parameters vary randomly in the range: M12 = 300 ~ 30006*61^, / = 400 ~ 3000GeV^, 
Xl = 0.1 ~ 0.7. We can see the plots of are entirely between the 2a and ?>a region 
of its experimental value. 

Finally, we discuss the Zhh couplings in the LHT model. In our calculation, we still 
consider the above three scenarios and let the parameters vary randomly in the range: 
M12 = 300 ~ 3000^6^, M3 = 300 ~ 3000Ger, / = 400 ~ 3000GeV", xl = 0.1 ~ 0.7, 
the numerical results are summarized in Figs. (7-8). We confirm the result of Ref. 18| . 
in which the correction from the mixing diagrams between t and to Zbb couplings 
is mainly on the (^1 and doesn't have the correct sign to alleviate the large deviation 
between theoretical predictions and experimental values. The plots scatter beyond the 3a 
region their experimental values are mainly caused by these couplings. Furthermore, the 
correction on the g^^ is very small. However, there is a little difference when we consider 
the contributions involve other new particles. At this time, we can see part of the plots 
scatter in the 3a internal region of their experimental values, where the deviation of 
can be alleviated. Unfortunately, the correction on the (7^ is still very small and the plots 
still scatter near the 3a region of their experimental values so that the large deviation 
between theoretical predictions and experimental values can't be explained. The similar 
results are found on the g\ and gy. 
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IV. CONCLUSIONS 



In this paper,we studied the one-loop contributions of the new particles to the Rb and 
Ap^ for three different scenarios in the framework of the LHT model. Prom the scatter 
plots of Rb versus xl, the precision measurement data of R^ can give strong constraint on 
the Xl- Considering this constraint, we can see Rb isn't sensitive to the mass of the first 
two generation mirror quarks. The relevant parameters are weakly constrained by the 
precision measurement data of A^p^. In the given parameters space, the large deviation 
of Ap^ can't be explained reasonably. Prom our study, the LHT model can provide the 
correction to the and have small part of the parameter space to alleviate the deviation 
between theoretical predictions and experimental values. But the LHT model can't pro- 
vide the large correction to the so that the large deviation between the SM prediction 
predictions and experimental values of the Zbb couplings can't be alleviated substantially. 
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Appendix A: The expression of the renormahzation vertex f ^ 7 



20| 
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Appendix B: The explicit expressions of the 5gj^^ 

They can be represented in form of 1-point, 2-point and 3-point standard functions 
A, Bq, Bi,Cij. Here and Pb are outgoing. In all expressions, the mass of b-quark 
is ignored. 

SQL = ^g'CliVHaUVHahml^CS 

+ - MlJ^Bo(-Pb,mai^,MA,)] 

11 1 d 



2 '^H^ 2 2 



1 d 1 1 

+ 2^d'^Kk - Mk)-g^M-Pb,ma^,Mz„)] - -S^[--ml^Boi-Pb,ma^,Mz„) 

- MlJ-^Bo{-Pb,m,^,MzJ]} 

1 d 1 1 

+ 2^«i,("^n|, ~ ^w„)^Poi-Pb,^yi^^Mw^)] - -S'^[--ml,^BQ{-Pb,m^^,Mw^) 



1 9' 
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+ACL - 2] + [^Bo{-Pb,ma^,MA,) + \{m\ - MlJ-^Bo{-Pb,ma^, Ma,) 
-ls^wBo(-Pb,m,i^,MA,) - Isliml^ - Ml^)-^Bo{-Pb,m,^,MA,)] - ^ + ^5^} 
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1 c) 
--Slr[mlBo{-Pb, mt, Mw,) + ml{ml - M^J-^Bo{-Pb, rut, Mw^)]} 

g'^x\ 12 
^j2(^C!^^ytb(ycKM)tbY^{~'^^wBi{-Pb, mt, Mwi) 

~S^w[Bo{-Pb.muMw^) + {ml - M^J^B^{-Pb,muMw^)]} 



+ 
4 
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